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i LNReDUCTLON 


A prediction interval is an interval that contains a 
future observation with a pre-specified probability. The 
limits of the interval are functions of known observations 
from a family of known distributions with given parameters. 
Though prediction intervals resemble confidence intervals, 
they differ from the latter conceptually. A confidence 
interval covers the value of a parameter of a distribution. 
A prediction interval, on the other hand, encloses the value 
of a random variable. 

A prediction interval for a single future observation 
and simultaneous intervals for independent and identically 
distributed samples from N(u,o7) were derived by Hahn [2]. 
For correlated samples with a prescribed correlation struc- 
ture, prediction intervals were developed by Choi [3]. It 
was shown by him that many of Hahn's prediction intervals 
are valid even for the case where the samples were correlated 
with the specified correlation structure. 

In this thesis, the observations are assumed to belong to 
a multivariate normal family with mean vector 


LS Gh sly DTD 


igo <il 


miGucOmeetat.on matrix 


V = 6 l + oo Cine mia) 
eee 50°) 1 








Let KX (t=...-2,-1,0,1,2...) be a stationary Gaussian 
markov process of order one with E(X,) = u for all t and 
the covariance between x and ars = 6-8 Pots lir teand k.. 
Stochastic processes of this type are used to model certain 
types of oceanographic and metereological phenomena and 
some aspects of the stock market. Any finite sub process 
will then belong to the above multivariate normal family. 
Chapter 2 deals with prediction intervals from a normal 
distribution and contains a summary of known results and 
a review of general linear hypothesis. In chapter 3 a 
prediction interval for a future observation in a first order 
markov process is derived and the properties of this interval 
are examined. The effect of various parameters such as 
variance, correlation and sample size on the prediction 
interval are also discussed in this chapter. Chapter 4 
contains a discussion of a simulation to generate serially 
correlated random variables, which are used to empirically 
verify the theoretical conclusions drawn. The results of 
the simulation are presented in chapter 5. Extensions of 


the results to a linear trend model are discussed in chapter 6. 





It. PREDICTION INTERVALS FOR NORMAL DISTRIBUTION 


A. SUMMARY OF KNOWN RESULTS 

Consider a random sample of size N=n;+n- drawn from a 
normal distribution N(u,o*). The samples n, form the group 
of initial samples and Ne form the group, called future 
samples. Based on the mean Xo and variance so of the 
initial samples, prediction Peels rots xX a sé were 
obtained by Hickman [6]. Hahn [2] derived prediction 
intervals for the mean and variance of the second sample 
and also simultaneous prediction intervals for the variance 
of each of k eee random samples of size Ne based on 
the information obtained from the initial sample. Hahn's 


[2] two sided 100 r % prediction interval to contain a 


single additional observation X is given by 
on ie 
Yeas te Pl he (Ole me) yan hues) Ss 


where Y is the sample mean, e- is the sample variance and 
fia), (ltr) /2}] is obtained from the tables of t distribution 
with (n-1) degrees of freedom. A 100 r % simultaneous 
prediction interval to contain X,, Xo, . XK, k future 
observations is given by 

Wart idk. 11,2) o 
where 

20 
R(k,n,x) = U(1+5)* 


The value of the function U can be obtained from the special 


tables given in reference 2. 





A prediction interval for a single future observation as 
well as simultaneous intervals for a specified number of 
future observations have been derived by Choi [3] for the 
case where the observations are correlated and belong to a 
multivariate normal distribution with mean 


ab 
m = (pps coc Ola 


mma cOVarlance matrix, 


V =¥(H+H') ta(1l-E) 
nxn Igbo sob eg! nxn nxn 
hy hy hy 
where oa ho hho ... hy 
nxn 
ho Dy by 
h,(i=l,2,...,n) and a@ are positive constants. E is an nxn 


matrix all of whose elements are unity. 
B. DEFINITIONS, NOTATIONS AND REVIEW OF SIMPLE LINEAR 
REGRESSION ANALYSIS 
A brief review of the results and notation in a simple 
linear regression model are presented below. These results 
are used in chapter 3 in deriving prediction intervals for 
a first order markov process. By way of notation, a 
capital letter with an underbar represents a vector or a 


matrix. Consider the simple linear regression model 


Y= X Bte (2-1) 
jobads ee x abel 


where eae N(O,0°L) 








mre an Sealar motation can be written as 
Yo ="4a + bx, a 5 = eee on. 


a: 


The least squares estimator for B is obtained by 


vib @ = (G2 > EO VCGS sea), 


that is B = Sey 
where = xtx 
‘3 jel 
and BeieNnCS, o 5 ). Ze) 


Beis the column vector Ca 


fee isewell known that the random variable 


¥ ~ N(XB, o°Z) 


rxé =i 
-l _ al 
-) ~ nse | 
-)?.X. nN 
cE 
where 
a D 
Sex = 2} (XX. = X) 
- a 
ie 


et So be the usual unbiased estimator of se given by 


eye b-S0 a = 2) , (2-3) 
S ay, 
mm@ese SYY = +» (Y. - Y) 








Let Y be the predicted value of the random variable at a 
specified value of the independent variable x. 
A 100(1-a)% prediction interval for the dependent vari- 


flowy eecan be Obtained from the relation 


P ___ ee <t]f?-= 1-a 2a) 
—, 2 
* 14 (&-x) 
Come: + = + SRY 


Equation (2-4) will be used in deriving a prediction interval 


bem x in chapter 3. 


C. STATIONARY GAUSSIAN MARKOV PROCESSES 

Let X, (tc = ...-2,-1,0,1,2,...) be a stationary gaussian 
markov process, 1.e., an autoregressive stochastic process, 
of order one with E(X,) =m and covariance between X and 

eek 
Pei, 9 - 

To derive the prediction interval for a single future 
observation Kony] based on the observations up to Koy» the 
results of Ogawara [1] will be used. He has shown that 
when the values of Xo, _7 (k=1,2,..,ntl) are fixed, Xoy 
(k=1,2,...,n) are mutually independent. The conditional 


probability densities are given by 





i 1 2 
meas | X x ) = EXP - ay fx : caren] 
mle ic 2+ 2 2k 
V 270. 20, 


5) 


where 


ane Oy) Clea 


0) 
ll 


oO 
ll 


‘ey Glega 
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oo = a8 (1np)/ (1t9*) 
X= (Koy Ft Xap yy) / 2 


Peprediction interval for Xonay will be derived in chapter 
3 based on this conditional distribution of the even numbered 


observations conditioned on the odd numbered ones. 


IEE 





Cine DLC LON UNteRvVAb FOR A FUTURE 
OBSERVATION IN SERIALLY CORRELATED NORMAL SAMPLES 

A prediction interval to contain a single future obser- 
vation for samples from a multivariate normal distribution 
with a correlation matrix of the form defined in equation 
(1-1), is derived in this chapter. In section A a conditional 
prediction interval to contain a single additional observation 
based on the conditional distribution of Xo, is obtained. 
Section B deals with the properties of the prediction interval. 
In section C the dependence of the prediction interval on 


sample size, 9 and sigma are discussed. 


A. DERIVATION OF THE PREDICTION INTERVAL 


Assume a sample of observations, X,, Xo, X, -» Xop ly? 
Xo, from a normal distribution with mean zero and covariance 
; Wes: 
matrix V such that Cov (X,, X,,,) = 9°. 


The conditional probability densities, as indicated in 


equation (2-5), are 


el ail _ sal 2 
E(Koe| Xoe-a> Roe) > a BP ey: ir catox.| 
AO ‘o1 
O O 
with conditional mean 
E (Xo, |X) = atbX) , 
Sr = (X54 Xo) /2 - (3-1) 


E(X5, |X) Gonistdenedsaqsed tunceton Of Ky is callled the 


regression function of Koy on xy. Graphically, it represents 


TZ 





the locus of the center of gravity of the conditional random 
e t . ' 
variables Xone | x as a function of x. 
As a result of Ogawara's Theorem the following conditional 


simple linear regression model may be assumed, 


at bxX/ te 


Koy = me Bp 


[eee es, TY (3°2) 


where Shem NCO, 07). 
The maximum likelihood estimators of the parameters a, b 


and oe are given by 





bo = oX1/SXX 
a= Xo, ~ b x 
7. Bev 
6 = (SYY - b SXX)/(n-2) 
where 
ah oe 1 X1 7X, . XatX 5 i A Kop _ pt Xop4y 
i in a Ds is 2 
Xp, = +(x, + x, + + x 
DK n a2 4 a: 2k : 
3 — 
SXX = D(X}? a= 
k=1 
S ' re | 
SXY = z x, Xo, - MX, Xop, > 
k=1 
n 
2 —- 2 
SHY = ( Xpx | - NM Xo, ; (3-3) 


At a known value of Xp the predictor x5, is given by 


nN 


Koy (3-4) 


= = ~ b xy 


Re 





A 100(1l-a)% prediction interval for the dependent variable 


Koy Ean e OpEaIneas trom the relation 


: =e Om 63° 5) 


From the above relation a prediction interval for Xoy,47 can 
be obtained noting that x, is a Mimeatashunerton Of Xopay- 


The inequality inside the brackets can be written as 


a eB : ie = + Oy - %) ; (3-6) 


Substituting for x,, = (xX. - b Xy) + b Xp and squaring both 


sides, the inequality (3-6) reduces to 


Sea Z te: 1 a iz 4 ee { = mee | 
(Xo, - Xo.) + BG - RYU - 2bCRop- Koy) CR, - XL) 


< to*(1 + 1) + to" (xy! : % 1) “/SXX 





, 272 z 
2 +29 See NG, = ae 
(5 - SKK (xp - x | - 2b(Xo, - Xo.) Cy - %) 


—- \2 Ze 1 
+ (Xo, - Xo) =~ to (i + <0 (3-7) 
This is a quadratic equation of the form 
’ =, 2 { A | 
A(X - xy) ~ B(x, - xy) +c s O. Co-3) 


14 





where 


ee te 
SEX 
B= -2b (Xo, - Xo) , 


i = 2 22 nr 
C = (Xo, = X.) - t oOo a (3-9) 


Thus, the probability statement in equation (3-5) is equiva- 


lent to 
PLACx, . mI) + B(x, - Xp) CaO = Lo =o % (3200) 


Pmeceaiction interval for x, Can now be obtained in terms of 


the roots of the quadratic equation (3°8) viz., 


-Biv B* - AAC 





(xy - x) — aera, Ye (Seq) 
the term under the radical sign simplifies to 
2 249 (xy, = Koy)” atl (52 "| 
B - 4AC = 4t 0 SXN — — b = SXX 


and hence 






D(x ; Ce 


Since x = (X5,_1 + ¥o,44)/2 » 


The unknown future value Xo] Cane pee we lecennin terms ot Ehe 
known Value as 


Kod] = D + E where 


1S) 








Z —; ~ i aoe), 
See x eb (Xo, - X>,)/(b - 








e9) 
li 


2*9 2-2 
‘ SS 32 Mle, ERO Le See 
2ta/ (Ko, ~ Xo.) /SXX + ——(b™ - mae) / 7 nie) 


(3°14) 
Thus a prediction interval for Koay 1s given by 


IL) > Boas Xora, $$ ID a> 1c; Se C375) 


where D and E are given by equations (3°13) and (3:14) respec- 
tively. The properties of this interval are discussed in 


Section B. 


B. PROPERTIES OF THE PREDICTION INTERVAL 
The prediction interval in the previous section was 


obtained by solving 
’ vl 2 ' v1 * 
A(X, = x) ae BCX, - x) eo 0 


Meee Enencocrinerents A, B and C are defined in (3-9). The 
roots of the equation and hence the form of the prediction 
interval will depend on the relative sizes of A, B and C. 

The effect ot these coefficients on the prediction 
interval will now be examined breaking in to four all in- 
clusive cases. 


Z 


Gase 1. A > 0 and B™ - GAC > 0 
A> O »> b aor 


16 





Or 7 


oo (Cy 
ESO 
boSKX _ 
Let Za = F (3216) 


Peme) GAC implies that 


2° 2 
4 52 (x = $3 y2 4 (62 aint i ( X 
2x 7~ *2,° > - Bx) So, 7 Xo) - 


EG : to Z 














ee Sy teem 
(<2 fs? | £25 ? 225? _ 
or b° - SX (n+1) > - her (Xop - Xo) 
Mee, | On 82k ~ X 2%? 
5p, LL ae 
Ee} n+1 = ay 
<2 (Xn, - X 2 
Ae DESO eee 2 2k 2k 
ag) n+l 
| Sa 2 (ce) 
To satisfy the conditions A > O and a > GAC it is necessary 
emat 
(X94, - Ky)? 
k 
Peeeieand fF > 1 — — an (3-16a) 


Since the quantity (Xo, 7 Ky.) */ too" is always positive the 
conditions in equations (3-l6a) are equivalent to the condition 


foie in this case the prediction interval for Koya] i a ae 


sey 








be of the form 


| (D - E) , (D + E) | 


This type of two-sided interval is referred to as TYPE l, 


where D and E are as defined in (3°13) and (3:14) respectively. 





ase 2. B- ——~oAGe> QO and A < 0 
., eo? 
ee) see Dox 
Or Da °< SXY 
pe 
b SKX 
Or ic) — 1 
eG 
free en < Os F< |. Cry) 
~ Z 
Gon =m A.) 
> So DUNG Ss 0) > 10 Sa ee eS 2S (3-18) 
n+1 eD 
8) 
2 


n  “%ox~ 7 2%) 


mnaess if F < 1 and F 21 - Posy] La 


the resulting interval will be of the forn, 


(- ©, D+E] [D - E, © ) (3-19) 
and will be called a TYPE 2 interval. 


Gaise Ss. If Be Pe NGO lanalsA + /© the interval is calied 


mers. 3. 


18 





The roots of the quadratic are complex and the prediction 
interval must be taken as (- ~, +»). This situation occurs 


when BA - GAC < 0O. 


— Z 
Z n (X547 Xo) 
Be - 4AC < O > F< 1 - ey a (3-20) 


Thus a type 3 interval results when 


ey) 
(OC ee tates) 
n vara 2k 
eo aT IES 
|e) 
PACE 4. AZ=0 
9 
= oo = (3.21) 
le Gs 


Thus, when F = 1, the interval is of type 4 and is of the 


norm (|G,> ) , 

x - x 
ntl) a 2k 2k (3022) 
a b 


where G= 2X). - Xo,_4 7 lo) SyOe 1 
To summarize the results, the quadratic equation that 


leads to the required prediction interval is of the form 
Meee) > + BC" = x') + C= 0 
k k Ks k 


and four different types of intervals can result depending 
on the value of 
bo SXX 
ee! 


Mae LOur cases are 


Ws 





i) when F > 1, the interval is of the form 
ii) when F < 1 and F 2 lL —7T —— the interval 
is of type 2 and is of the form (- ~, DtE],[D-E,~) 


— 2 
(X5, = X51) 
fet) when F < 1-= a , the interval is of 


BE 0) 
type 3 and is of the form (-~,+-) 


iv) when F = 1, the interval is of type 4 i.e., [G,~). 


At a given level of significance a, from the relation 
(3-5) and the conclusions drawn in this section, the following 
relation should hold: 

Ly 

__ [occurrence of type i interval] 


otha rath interval contains the predicted value] 


=l-a CBnez3)) 


The above conclusions are verified with simulated 
samples in chapter 5. 
Section C examines the distribution of the quantity F 


defined in equation (3.16). 
See OP oPREBULLION OF ‘F' 


The occurrence of any particular type of interval depends 
Oem the value 
b¢ § 
FS yy 
ELS 
The distribution of F can be shown to be a noncentral F 


feseripution as follows. 


20 





Mmeem (2°2), it can be shown that is 1s N(b,o7/SXxX) 


(b-b) 
or ———— rv NCO,1) , 
o/¥SXX 
a 2 
or &-b) Nv x7 (1), 
oo Ok 
b-- SXX 2 
and i ae ay x (1, r) : G24) 
fe) 
b? SxXX 
where A = 5 ; 


is the noncentrality parameter. 


fegom (2-3) the distribution of 


a? 
{nit}o is \2 (n-2) . (3-25) 


Oo 
The ratio of (3-24) and (3-25) devided by their respective 
degress of freedom results in a noncentral F [5] statistic 


with degrees of freedom (1, n-2) and noncentrality parameter 


me Chat is 
b*. SXX 2 
0) = t FF. F(1, n-2, A) , C3 6) 
Oo 


and the noncentrality parameter 


pax 
a 


Oo 


r 


Substituting for b from equation (2:5), the noncentrality 


parameter X can be written as 


ZX 





2 Lo Sp 0,4 
i = Se (3-27) 
(1+p*) gf 


D. NUMERICAL EXAMPLE 
The Model developed in section A is tested on the data 

obtained from Dow-Jones monthly averages to predict a future 
monthly average, based on the averages of the previous months. 
In most of the cases, the future value was contained in the 
prediction interval. To illustrate an example the data for 
the years 1966 and 67 are given below. The prediction 
intervals are computed with the help of the APL program 


PREDICT, whose listing is given on page 68. 


1966 i ey 
Month Ended Dee J Ind Month Ended D = J ind 
January 31 OS so January 31 849 .87 
February 28 951.89 February 28 35) 5) 5 Sai 
March sul 924.77 March Sil 865.98 
April 29 SN 8h 815 (exe! April 28 897.05 
May 31 884.07 May S580 eh 2k ONS) 
June 30 S7.0.00 June 30 860.26 
July 29 847.38 july. 31 904.24 
August Syl 788.41 August Sal SMO 28. 
Sept 30 WHE V2 Sept 29 926.66 
October 31 oC TE Ol OcEober jo! 879.74 
November 30 SNL BS, November 30 875.81 
December 30 Tes ayeeM) December 29 SHOVE al 


V7 





The above data are stored in the variable DMA. The pro- 
gram Predict is run with different sample sizes and the out- 
put is shown on pages 24 and 25. 

Using 14 monthly averages from January '66 to February 
'67, prediction intervals tor March '67 are computed and the 
exact value of 865.98 is contained in the interval. 

Using given values up to April ‘67, the prediction in- 
tervals for May ‘6/7, also contain the true value. With a 
sample of size 18, the prediction intervals for the nine- 
teenth month i.e., July '6/7 are computed and they contain the 
true value. Though the average of September ‘67 falls outside 
the interval, the average of November ‘67 is contained within 
the interval. 

The reduction in the length of prediction interval with 
increase in sample size, is supported by the above example. 

The data are tested for the required correlation structure 
and found to fulfill the requirements partially. The statistics 


perenew Gata are 


Mean = $870.83 
Standard deviation = oO) 
Gomer COCEEILCIENE /— = 0.82 


Ze 





PREDICT 14tDMA 


ie fF STATISTIC [5 

THE LOWER PREDICTION LIMIT 
iis UPPER PREDICTION LIMIT 
PP GeNGlH OF THE INTERVAL 
TRUE VALUE OF THE AVERAGE 


PREDICT 16+DMA 


mit SlATESTIC IS 
EitwnOwek PREDICTION LIMIT 
item PREDICTION LIMIT 
THE LENGTH OF THE INTERVAL 
TRUE VALUE OF THE AVERAGE 


PREDLC) £364 DMA 


Gites tAliSrTiCc Is 

he eeOWwiR PREDECTION LIMIT 
fimevPPBR PREDICTION LIMIT 
THE LENGTH OF THE INTERVAL 
TRUE VALUE OF THE AVERAGE 


PREDICT 20+DMA 


PibetaolArISTIC [5 

THE LOWER PREDICTION LIMIT 
THE UPPER PREDICTION LIMIT 
THE LENGTH OF THE INTERVAL 
TRUE VALUE OF THE AVERAGE 


24 


a, 1320/3678 
J ie VEZ SO 28 
1080,631201 
S155), Skis) Si iLtaits: 
865.98 


5.489038103 
SSS) 5 SSIS) Sioa’, 
900. 5564036 
S02. U6Z28487 
ele) 4ts DIS, 


Det SO ot 243 
THOtok OS 

J GIOIS) , 2b ONILeIENe 
B00hs 979965 
904.24 


ee no Soa 
Sif) 8, SNA SMas Pals 
864.8705316 
291.277549 
7210). elle 





PREDICT 22+DMA 


mee F SlATIsiIic IS = 6.994239226 
Pee LOWER PREDECTION LIMIT = Sees ILI eles 
ines e ER PREDICTION LIMIT = 899. 7280383 
THE LENGTH OF THE INTERVAL = ZOOmole oi 
TRUE VALUE OF THE AVERAGE = Via) pie te 
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IV= SIMULATION OF SERIALLY CORRELATED SAMPLES 


To study the effect of op, o and n on the type of pre- 
diction interval obtained it is necessary to generate random 
variables with the required variance covariance matrix. 


These can be generated as linear combinations of standard 


normal random variates. 


A. SIMULATION FROM STANDARD NORMAL SAMPLES 


Suppose Z4: Zorsees Z. are random samples from N(O, 1). 
To get serially correlated random variates Xy, Xy,-+-,K,, 
let Ky = 24 
Xo = C,Z,+C,Z,, 
C2, +C5Z5+ we +012, : (4-1) 
The coefficients Cy. -. 4,0, are then determined by solving 


equations of the form, 


esd 
Cov (X, ae) =o for allt and k from’ 1 to mn. 


Consider the first equation 


x, = 2 


1 1 


E(Xy) = E(Z,) = 0 


Var Ca Var (Z,) = 1 Cae?) 
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From the second equation 


E(X5) = E(C,2,+C,2Z,) 


C,E(Z,)+C, E(Z,) = 0 (4-3) 
since Z4 and Z5 are independent; 


Var (Xy) = Var (C,Z,+ C5Z5) 


= = Var Z, + cs Var Zo, i.e., 
B+c,* = 1 (4-4) 
Cov(X,, X,) = E [(X, - E(X,)) (XK, - £(K,))] 
= E [X,° Xp] 


= £E [Z, ° (C1 2,+C,2Z,) J 
= @§ fee ew) 
egal! Dag +e, 
Ca (Z 246 EZ. 2) = 
eal aes?) 
2 2 
Since Z, v N(O, 1), 2," v x" (1); 


also since Z4 and Z5 are independent, the above equation can 


be simplified to 


Cy =o C25) 
From 4.4 by substitution we get C, = lL - ae 
or =X, = pZ, +¥1-p* Z (4-6) 
Z il! 2 
or es (Ox, + 1 -o Zo (4:7) 








Simiiarly it can be shown that 
Xy = 0X, +y¥l-p* Z (4-8) 
5) ay 3 


The result can be generalised in the form 


X, = OX. 4 ty l - 0 Zs (4-9) 
ily = 4 ne 7 a 
ae = lt + a" ea eee yi - 2 7 (4-10) 
i P te p mn CE a 3 
ee ee ee. | TY 


The results obtained in (4:10) can be expressed in the follow- 


PPeemachEx notation: 


X i i oe (oN) 


nxl nxn nxl 


where i SO Xa, TOK ie ‘ 
— L Z n 
msc 1 
= le 
Z aa (25, Zo; ry s s 9 Z) r 
nxl 


and A is a nxn matrix of the form 


i 0 Chee. Sere 0 


Gari) 
ss o¢ o¥l - 0 yl - a eg. 


It a random variable Z is distributed normally with 
nxl 
mean yp and covariance V and if A is an nxn matrix of rank n, 


then the vector 


28 






is distributed as multivariate normal with mean vector Au 
and covariance A V Ar, As a result, the random variable X 


(4-11) is distributed as multivariate normal with the re- 


quired covariance matrix V. Here 
Z~vN(O, I) , 


where I is an nxn identity matrix. 


Hence, Xv N(O, AL a‘) 


We now show that AI al = A L 


PSeasiadeci xX wor Ehe required 


mommeccor n = 4) 


iL 0 0 0 
fm = O /\-9* 0 0 
4x4 
7 a) ee 0 


3 ‘ l-o p¥l-o /1-9" 


Z 3 
iL 0 0 0 
al. | 9 Yi-p oY¥i-p 0-0" 
x4 
pe 0) 0 Y1-o pv¥l-o 


0 0 0 ee 
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oa hace) 








2 2 
pepe 1 Pelee) pee aes) 


3 
p? e MiGiag 5 peta leo+ ee Clseo) ots: Clee ea (Clore 


9° meter (le) Bote" (t-p-)4o(1-0°) j FUL) Coeuwal aly 


which reduces to 


L 0 ie ae 
p 1 p a 
— 
AA 9 
fe 0 i fe 
a ~~ oe 


The above result can be generalized for any value of n. 


B. TESTING FOR NORMALITY AND CORRELATION 
The simulated samples are tested for normality and the 
correlation structure. The simulated samples are tested by 
standard Kolmogorov-Smirnov test. The results showed that 
samples do appear to be from the specified normal distribution. 
The correlation between Xx. and X, is given by 


COV OC a) 
R= . eee 13) 


Vilibe: Gaunt iebcn.. 
i J 


Confidence limits for correlation coefficient p are developed 


using the approximation. 


oy) 








Ree te l+o 1 
Z, 5 in ik N(% ln loaf eel, 


A 95% confidence interval for p is given by 
(698-1) / (67841) < 9 < (e2P-1)/(e7P 41), (4.14) 


where 


ro 
\O 
OV 
ro 
\O 
OV 





and b = Z 


+ 





al 
Lo 
2 
i 
Lo 


The simulated samples are also tested for the required 
correlation structure by examining the confidence limits for 
correlation between (X,, mY using equation (4-14). The 
results showed that the samples do appear to have the pre- 


Serapea correlation structure. 
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V. RESULTS OF SIMULATION 


In this chapter, using simulated samples, the probabili- 
ties of occurrence of various types of intervals and their 
dependence on the value of 0, n and o are studied. The pro- 
gram PROBS FORTRAN, whose listing is given on pages 63-66 
computes the probabilities for the four types of intervals. 
For given values of 9, n and o the program generates one 
thousand samples and records the frequency of occurrence of 
each type of interval. In each case the program tests whether 
the predicted value is contained within the interval. Finally, 
for a given level of significance a the relation (3°23) is 
verified. 

Choosing a level of significance a = 0,05, the simulation 
was run for values of o ranging from 0.1 to 0.9, for sample 
mezec trom 10 to 50 and for sigma varying from 1 to 5. For 
ao = 0.1 and 0.02 the program is executed to study the effect 
of the level of significance on the results. 

The computer output is shown on pages 43 to 57, Tables 
I to VI show the variation of probabilities of types of 
intervals with change in a at different sample sizes, keeping 
ao and o constant. Tables VII to IX show the variation of 
probabilities of types of intervals with a change in 0 at 
values of a = U.1. Tables X to XII show the variation of 
probabilities of types of intervals with a change in 09, 


keeping a = 0.02. 
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im=eaples ALTil to XV, the variation of the probabilities 
of types of intervals are shown at a standard deviation (oc) 
of 3.0, whereas in Tables XVI to XVIII, the value of standard 
deviation is altered to 5.0. Tables XIX to XXVII depict the 
variation of the types of intervals with change in sample 
sizes at different values of p and o, keeping a at 0.05. 
Finally, Tables XXVIII to XXX show the variation of proba- 
bility of types of intervals with change in sigma at different 
values of o and n. 

In all these tables the first column indicates the value 
of the parameter. Columns 2, 3 and 4 give the probabilities 
of occurrence of types of intervals. Columns 5 and 6 show the 
probability that the predicted value is contained in the 
interval. Since a type 3 interval is of the form (-~,+~) 
the probability of Xop41 £0 be contained in the interval is 
ome finally columm 7 verifies the relation (3.23). The 
attached graphs on pages 58 to 62 show the trend in the pro- 
bability of occurence of the three types of intervals with 
changes in parameters like oe, o, n and a. 

The outcome of the simulation confirmed the following 
theoretical conclusions: 

i) The probability of type 1 interval increases with an 
increase in correlation coefficient o, all other parameters 
remaining constant. 

ii) The probability of type 1 interval increases with an 


increase in sample size, keeping all other parameters constant. 
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iii) The probabilities of types 2 and 3 decrease with 
an increase in correlation coefficient or sample size. 

iv) The probabilities for the four types of intervals 
do not change significantly with change in sigma. 

v) The probability of a type 1 interval increases with 
an increase in level of significance a, all other parameters 
remaining constant. 

Section 6 deals with the application to linear trend 


model. 
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VI. APPLICATION TO LINEAR TREND MODEL 


In this chapter, the results of chapter 3 are extended 
to stochastic processes that follow a linear trend model. 


Let 1S be a discrete stochastic process satisfying 


ae i (6-1) 


where = a + Bn , (6°2) 


and En is a stationary process of normal variates that 


satisfy the following: 


Bie} = 0 
Le 
varie} =O 
Tek 
Cor{e_,€.44)} =o 0 Oe (6-3) 


Such a process is called a model for linear trend. Krishnaih 
and Murthy [7] showed that for the above model the conditional 
distributions of X,, given X5,_4 and X5,,), are independent 


mommal with 


E (Xo, | Xop-y> Xonezl = 81 * BoK + 83%, (6-4) 

d Var (X., | X SS oO poy eS (6-5) 
cute ee 2k) 2k-1* ~2k+1 
where » eS n,n Xo 4) /2 
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imeseclion A, a prediction interval for Kova is 


obtained knowing the values up to Koy for the above model. 


ee fr aeDLCTION INTERVAL FOR Xovay 
Using regression analysis, a prediction interval for 
Koay can be obtained tor the Linear trend model, knowing 


the values up to Koy. Equation (6-4) can be written as 


Koy = 84 a Bok sta BX) + EL: r= 1) eee Te 
iain macrix notation, 
Bo. - X& B+ « . €¥NO, OD . (6:6) 
max | nx3 3x1 no 


The design matrix X is of the form: 


mS 
1 1 (X, +X) /2 
aL 2 (X3+K.) /2 
1 3 (Xo+X2) /2 
Lon Ky Rong) /? 


and B is the column vector (84. Bo 83) 
3x1 


The estimator B can be obtained from the normal equations 
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xe xy pe xt x . 
5 X ok 


3xn nx3 3xl 3xn nxl 


Mime beE 
= BSS XX, > 


n ag! De 
n 


where 3 X° X = Het Ea Enx! (O27) 
Speigl i9y:<8 3 
ie} 0D,G ee Ye 
n n n 


Since S can be shown to be nonsingular and of rank 3, Sas 


exists and 


Jip it 
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Bu N (B, ee (6-8) 


For a given value of X', = [1kx 
x 


the predicted value X = xX’ 
bg a 


a t Zz t == Ll ' a 
and Xo, » N{X, B, o (X08 “(Gy D7} 


2k 


It can be shown that the residual 


- Lape D | 


where 


SR 


S nig = 
lied m3 axeor OxXk 
As an example, let X = (X}, X59, ---»%4Q) 


Here K = 5 andn = 4. 


The vector X45, = (9, Xy, Xe: Xe) 
4x1 


5). 








1 iL (x, +x3)/2 
iE Z (xXgtx_) /2 
and the design matrix X = 
4x3 i 3 Gere) 2 
5 
1 4G (X54KQ) /2 


Ay 


=> iain S (xXgtx,,)/2] 


The value of X71] Ponce wh Ec pLedi eChlonainterval for X14 


can be developed as shown below. 


From equation (6-9) 


an 


x - x 
2k 2@k yO CO, 1) as) 
Cy = 
xe - x R 
and 2k 2k pO Fe (6-11) 


Cee /a-3 


The quantity t has a student's t distribution with (n-3) 
degrees of freedom, where RS is the usual unbiased estimator 
ong 5°. 


At a given level of significance a , 


Gr - . ) vn-3 
ct) = b= © 6-22) 


Ro Ho) Se 5° 


The term in the brackets of the left hand side of equa- 


tion (6:12) can be written as, 
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a PD fy 
and 


Les s1+(ikxti(s tyre xy 


The expression on the right hand side can be simplified to 


Ene. £Orm 


ft f Z 
ee ik OX 


Mine coefficients eon a and Yy can be computed since k and 
g~t are known. 
Squaring both sides of equation (6-13) and substituting for 
Xo, and 5 it can be simplified to, 
2 
f 2 a/2 O ? m a 
<Xp BTR gk-B 4%) S$ aa VG ATV OX ? 


or 


2 RZ -¢ R-y 
aD, Fal 2e O ' _ = BLD Ora ' 
(8 Ae ew mi + (2 (XB -Byk) Bg- —ALE_2 *) xy 
2 R-y 
ean. a Fa/ 2 oo 
The above inequality reduces to 
ig ! 6-14 
Peo Bx + C x 0 ; ( ) 
where 
ae. 
A a B 3 por Yot ar Ré/(n- 3) , 
2 2 
nee A te my. 
Z pee cot 
Beppe 69) 3 5 


a, 











2 Zz 


es en 2 Leap ole 

c — (X5p7 By 85k) a n- (6 ‘ ILS) 

The quadratic equation (6°14) can be solved for Xo Mist ric 
ene rormula. 

xi = 2BEPB*-4AC ge gy 

k 2A 
Simce xy = 2(Xo,_4+%o,44)> Ehe predictlonsinterval for 
Xop.4, Can be written as 

D-E < Xo] < DteE Ca say), 
where 

Da B/A 

E = /(B2-GAC)/A CGral7) 


Note: If A=0, then no interval exists; also when 34 < GAC 
E becomes imaginary and the interval is taken to be (-~,+-)., 
In section A, a numerical example is given to verify the 


model developed here, with the help of simulated samples. 


Eee NW RICAL EXAMPLE 

In order to compute the prediction intervals for a future 
value aoe samples satisfying a linear trend model an APL 
program LINTREND has been developed. The listing of this 
program is given on page 6/7. Error terms having the correla- 
tion structure (6:3) are given as input and for given values 


of a and 8, the program develops the vectors Koy and X 
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a ee LF ee OTe Oe 


For a sample of size 18, the output of the above program is 


given below. 


X le e004, S507, 47000, 327595. 5.984. Oe, 
a EVE 8 cekhe 
X ae | 2. J ee ene a O13) 2.685, 50570, 7.84909. 017. 
aa 8.416)+ 

1 sh 2.094 

ih 2 Se 

1 3 4.056 

Miesdesign matrix X = /|1 4 5/58. 
exs 

l 5 5.984 

1 6 i. vas 

iL y ioe 

1 8 8.388 


For the model (6-4), 


B=(-0.147, -0.597, 1.566) 
eae! 


The coefficients of the quadratic equation (6-14) 


2 


| 
ro) 


A(X, ) ar B(X) ae eS 


are given by 
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A = 1.443 
B = 35.801 
Ge=)- 72-382 


and the values of D and E are 
Dee 15.955 
E = 28.561 
poem the equation (6-17) , 
the lower limit = - 12.606 
the upper limit = 44,516 
The true value of Xi g ClO 515) tous Comeained in tne 


mmeerval. 


42 








I. VARIA™ION OF TYPES CF INTERVALS WITH ROW 


Severe SIZE 14 
SIGMA = 1.0 
meet 9 SIGNIFICANCE = 0.05 


— ee ee ee ee ee eee ee eee ee ee ee ee sp en ee eee ee ee eee ee eee ee a ee 


ROW Mueeis TYPE 2 TYPE 3° PREE OL. PeRGa 62 “CON REC 


= ee SS ce ee es ee es eee es ee es wee ee eee es ee eee ae pe ee ee ee eee eee eee eS ee 


= ee eee ee ee ce ee ee ee es eee ee ee ee ee ss ee ee eee ee ee eee ee ee 


owen NOP TYPES OF INTERVAES WITH &2OW 


Seeeele Size = 22 
Sos = 1.vV 
Mevieie cr SloNiFICANCE = 0.05 


ROW mare) Y PG 2 YO. JPR Oe doe CON REG 
GO eLUU C2452 0. 048 Co. 46C OF SS) a) C2438 die oe 
Ce 30C C.516 O. Cl? 0.407 0.924% 0.714 Ore oi7 
Je 500 0.483 0.110 0.407 0.909 9.318 0.936 
0-790 O.791 0.071 OPA. ite. Oe oa4 0.873 OR me 
0-900 C.902 0.050 0.048 QO. S47 Ce 720 0. S42 


—— << eee eee eee ee ee ee ee ee ee ee ee ee ee eee ee oe = ee eee eee ee eee ee ee eee ee eee 


CON REG = CONFIDENCE REGION 


4.3 








Ili. VARIA™ION SF TYPES CF INTERVALS WITH ROW 
Sartre ke SIZE =7>0 


SIGMA =e 
LEVEL OF SIGNIFICANCE = 0.05 


ee ee ee ee ee ee ee ee oe SSP oe eee a eee ee = oe a ee ee eee ee ee es ee ee eee ee ee ee 


IVe VARIATION OF TYPES OF INTERVALS WITH ROW 
SirrreLe SIZE = 34% 


SIGMA = 1.0 
Peeve sSioNIFICANCE = 0.05 


> ee ee ee ee ee ee ce ese ee ee ee eee 


ROW ieee | TVFE 2 j(WVRE 3 PROS | “PRM 2°CON PEG 


—_—> = sep se ewe oe cep cee oe ce es swe es se em ee eee Se eee eee 


a << <p ue cote cole es ee es ee ee es ee es se i ee se se ee ee eo es eee ee ee A ee a ese me eee 
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Ve VARIATION 9*= TYPES OF INTERVALS WITH ROW 


oer ee SIZE = 38 

SIGMA = 1.0 

Bevel OF SIGNIFICANCE = 0.05 

RUW ares t  “YFE 2) *YPE 3 ~ PROG 1” PROB 2 CON- REG 
QelD0 C.703 0.020 0.267 0. S4S C.467 0. S46 
0.30C C. 809 O. C36 OA ilesye) 0.95% Os205 Cae 6 
0.500 0.847 0.041 O aie2 0.956 0.854 Onsa7 
Je f00 0.966 OrOL > 0.019 Oe Sys) O27s5 Pie on 
0-900 C.994 O.C0% 0-002 0.957 1.000 OSs fi 

VI. VAR TATION OF TYPES OF INTERVALS WITH ROW 

Be PLE SIZE = +2 

SIGMA = 1.0 

LEVEL CF SIGNIFICANCE = 0.05 

ROW ace y G& 2 TYPE Sees reese) 2 CON [EC 


oe ee ce ee ce ce ee eS eee es ee se ce ee em ee eee ee SS Ss eS Se 


= <b ce co cee cee cee cee eee eee es ces ee es ee ee ee ee ee es ee ee ee ees es ee es ee ee ee ee ee ee ee eee ee 
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Vil. VARIATION OF TYPES OF INTERVALS WITH ROW 


perrle SIZE L4% 
SIGMA = 1.0 
Bevel €F SIGNIFICANCE = 0.10 


ea ee ee ee ee ee ee ee ee ee ee ee eee es ei ee ee eo eee ee ee ee > <> ce ee eee ee eee ee ee ee eee oe 


ROW ier Lo TY ee 2. TYPE 3 PROB 1 “PROB “2°CON REG 
Ce10C Go2zo9 QO. 111 Cie. € QO. 879 0.468 0.506 
02300 0.441 0-108 0.451 0.846 0.567 0.896 
0.500 C380 Orso 0.482 0.858 Osi 2 OseeGs 
0.700 0-501 0.138 C.361 O. 858 Camo 0 C2c0C 
C.90C Ce. 666 O- 114 Ce22C 0. 883 0.868 0.507 


VITIT VARIATION O- TYPES CF INTERVALS WITH RCW 
Seve Le SIZE 22 

SIGMA = 1.0 

mee, ©F SIGNIFICANCE = 0.10 


le <P eee eee ee ee ee eee ee ee ee eee eee ee ee ee oe eee eee ee ee ee 


ROW levee 1YPE 2 TY@E 3 PRO@B Lt PRESB 2 COUN REG 
0-100 0.524 O.C77 06399 0.876 0.468 0.&94 
0.300 Osa 9 O 01 0.300 0.873 0. 564% Gorse 
0.500 0.6517 Oem 2 0.261 QO. 861 Ce. 154 0. €84% 
QO. 70C C.853 O. C66 0.081 0.885 Oe ee) 0.837 
0.900 0.941 0. 023 O.Gs1 0.891 0.786 0.891 


> <p ee ee ee ee ee ee ee es ee ee ee ee ee es eee ee ee eee ce ee eee ee ee ee ee ee ee ee eee ee =e ee 











IXe VARIATION @ TYPES CF INTERVALS WITH 2CW 


Serer = SIZE Bh 
STGMA = aoe) 
RevGle 2F SIGNIFICANCE = 0.10 


ee es ee ee ee ee ee ee ee es ee ee ee ee ee ee ee es ee ee ee ee ee ee ee ee ee ee ee es ee ee =e <ee 


RIW Wiel 1¥ee 2° TYPE 3 PROB 1  FRGB 2 :°CEN REC 


Pe eee ee ee ee ee ee ee eee ee cee ee ee ee es a ee ee eee ee eee ee ee ee ee ee ee 


> ab <P ae oe ee Pe ee es ee ee ee eee eee ee ee <P ee > a a ee em ede ee —— amp ae ab ou 


Memvanioata EN OF TYPES OF INTERVAES WITH & OW 


Seemeke SIZE = 14 
SIGMA = 1.0 
meee OF SIGNIFICANTE = 0.02 


ROW Tveweel. TYPE 2 TYPE 3 PREB L- FRa6b 2 CON REG 


_—- << ae oe oo <ee we e  es e ee ee  es e  e e ee ee ee  e ap <p <b ep Dc cc ee ce ce ee eee 
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Xie VARIATICN CF 7YSES OF INTERVALS WITH 20W 


Semele SIZE ae 
SIGMA = 1.0 
meee DF SIGNIFICANZE = 0.02 


—_ ae <p om ee SS ee es ee es es ee eee eee ee ee ee ee ee ee ee ee ee es eee ee ee ee ee ee ee ee ee eo ee ee ee ee eee oe 


XIi.e VARIATION OF TYPES DEF INTERVALS WITH ROW 


peer ie SIZE 70) 
> PEMA = 1.0 
Rewicee OF SIGNIFICANCE = 0.02 


RQW ape — Yre 2g TYPr 3 —~PRgs £ -PRGe 2 CON RES 


—2_ ae oe oe oe ce es ee ee ee eee ee ee eee ee ee es eee ee eee ee ae Se 2 Se ee A se eo es ee ee eee ee ee ee ee ee 
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XITIIWVARIATIIN OF TYPES OF INTERVALS WITH 20W 


Semreet SIZE = 
Sees = oe 0 


iw 


2 


Beyer OF SIGNIFICANCE = 0.05 

RCW (ogee OYFG 2 FYFE 3 PROS L PROB 2 CON REG 
02100 0.067 O.C/78 ORs ys) 0. 851 C. 500 Cao 
Ce 30C Cornig G22 5 Cro 6 Gs.o> 5 Ge/04 Oreo 
0.500 0.445 Groping 0.438 0.3899 0.829 O2355 
0.700 0. 736 0.083 OsLtS 0.938 arose C. S45 
0.90C C.SOl 0.051 0.048 O. 548 C. 941 Ca55 6 

Mivemevon ta 1iON OF FYPES OF INTERVALS WITH ROb 


ore le SIZE 


SIGMA 
Meese EF SIGNIFICANCE = 0.05 


Teee 1 


= 2 


= 


= 


Wy ee 


6 
- 0 


ENE Ge am, 


PROB oi 


Prope 2 COMNVREG 


<> ee ee ce ee ee ee oe ee ce ee ee ee ee ee ee ei ce ee ce ee ce ee ee ee ee ee ee ee ee 2 a oe 


Ge. 10C 
0.300 
0-500 
0./0C 
0.90C 


Ce 652 
Qa255 
Gi Asien! 
C. 841 
65 She 


O-. C43 
0.106 
0.110 
0.062 
O. C16 


OPE Reb 
0.539 
OS BAe 
0.097 
Ge022 
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OS o55 
OR C44 
0.9338 
0. S41i 
O. 946 


0.651 
0.698 
C3555 
Ge 50s 
Ueo30 


0.943 
C.S550 
C. C4S 
Cat 
0.947 








XV eVARTATION OF: TYPES OF INTERVALS WITH ROW 


Srowr ue SIZE 39 
SIGMA = 3.0 
mevet CF SIGNIFICANCE = 0.05 


a > ee gc eee eee“ Sc eee > <eee ce ee ee ee es ee ee ee ee ee ——— <P eee cee oe ee ee eee ee eee ee ee ee ee ee ee ee ee ee 


ROW TYPE 1 TYPE 2 TY®E 3 PROB 1 PRGB 2 CON REG 
C.10C G.075 C. C71 C.854 0.880 0.620 0.554 
0-300 0.374 O.1ll 0.5 Jeg22 0.721 0.S40 
0.500 Oso 74 GrsO3 1 0.235 0.942 C.857 Co. S4E 
0.70C C. &S7 0.034 C.06S 0.950 Ge 4 Oso 2 
Q.90C C2979 0. CO7 C.O014 0.952 1.000 Orso 


—— <> ie EP cee ie ee ee ee ee se se ee eee ieee ee ee ee es i i ee i i ee ee ee ee ee ee 


XVI .VARIATION OF TYPES CF INTERVALS WITH 20 


mameLleE SIZE = 22 

SIGMA = 3.0 

Meyveee fF SIGBIFICANCE = 0.05 

ROW ioe ivec 2. |YPE SPs eines CON REG 


a ee ee se ee ee ee es ee eee ce ee eee ee es es es es i Sc ec Ns SS 


> << aae> cop Oe ce ee ee ce ees ee ee ee es ee ee ee ee ee es ee ce ce es es ee ee eee ee es ee ee ee ee ee ee ee ee ee ee ee ee eee eee eee 
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XVII .VARIATION OF TYPES GF INTERVALS WITH 80W 
Serre LE SIZE = 26 

SIGMA = 3.0 

meee EF SIGNIFICANCE = 0.05 





ROW ree o 1 TPE SZ TYPE 3  PRGS 1 PROB 2 CON REG 
Je 100 a lek: 0.097 Q./10 02896 0.619 0.$43 
0-300 0255 Oey OAreRSe: 0.929 Ce 701 Ce =o 
Q.50C C2560 Oe lil OR BZ Gh eee Fs lei) 0.549 
Q.70C C. 841 OseGEZ 0.097 0.941 0.993 0.9544 
0.900 0.962 OL01s 0-022 0.947 0.933 0.548 


Peve VARIATION DF TYPES CF INTERVALS With ROW 
Somece SIZE 30 
SIGMA = 5.0 
ewe JF SIGNIFICANCE 


= 0.05 


2 <p oe ee ce SE ce ce ee ee ee es ei ee ee ee ee ee ae ee eee ee —b> a a we ee ee oe ee me eee ee 


a ee ee ee ee es ee ee ee es ee ee ee ee ee ee eee ee <> <> a ce ee ee com ee ce com oe ce oe ee eae ee ee ee 


Si 








XIXsVARIA-ION OF TYPES CF INTERVALS WITH N 
Caine. C OEFF . = 0.300 


SIGMA =o 
mepeGk OF SIGNIFICANCE = 0.05 


2 ame aap ep eee ee ee ee ee eee ee ep ee ee ee ee ee cee ee ee ee [eee «ee eee ae oD > oe ee ee ee eee ee ee ae ee eee ee 


N logee 1. TYPE 2 TYPE 3 PRB 1 FReB 2 CCN REG 
4 G 300 0.040 0.660 0. 867 C.-600 0. 44 

E Cao 42 Q.C80 OP A) Oe oS Cott) Cleese 

8 55 Shel i Gone on Cs256 0. 942 0.737 ORS) Ss 
1iAe, Gea le 0.077 0.407 0.924% 0.714 0 aS, 
LX Ge2o9 0.106 0.635 0.946 C.698 « $54 

14 C.861 O0.C33 0. 106 OQ. S40 Coal os C.54C 


XXe VAR TATION OF TYPES OF INTERVALS WITH N 


Cm. COEFF. = se © 
SIGMA = 1.0 
BeweleyCE SIGNIFICANCE = 0.05 


——p comp cule < ceee e cee cee e e ce es s e e es e ce e ee  ee ce e e e e  e ee eeee e e  e e e e e e s e e ee es e e ee 


N ieee), TYPE 2 TYPE 3S PRegery il Peas 2 CCN REG 


—_— a a ce ee ee ce ee ee es sees ee eee ces es ee ee ee ee ee ee ce ec ee ce ce es ee ee ee ee ee ee ee ee ce ee es 


LO C2483 O.110 C.407 
ke Ce. 561 
14 0.874 


D2 





XXI «VARIATION OF TYPES CF INTERVALS hiTH N 


Caren «<CUEFF. = 0.700 
SIGMA = 1.0 
Meee GF SIGNIFICANCE = 0.05 


> <p oe ee ee ee es ee ee ee ee ce ce ee al > PP ee ee ee ee eee ee ee ee ee ee eee eee eee ee ee ee 


N TYPE 1 TYPE 2 TYPE 3 PROB 1 PROS 2 CON REG 
0.308 Ole (OTe 0.620 0.396 0.861 OAS le 
6 0.359 OP ig Oe) 074 Oe 20 C.874 C. 54S 
C2586 QO. 095 C.315 0. S20 C. 884 C. S42 
10 Carol QO. C71 C.138 0.934 Osi tls OR as: 
L2 0-840 0.965 C2045 0.940 0.908 0.544 
14 C.91% 0.028 0.058 Oe >n Cees eee eS 


oe ee es ee es eee es es es ee ee ee ee ee ce ee es ee ee ce ee ee ee es ee ee ee es es ee es es ee cee ee ee oe 


Mette V AR TATION OF TYPES DF INTERVALS WITH N 


Gaivwe. CLEFF. =) 53010 
SIGMA = 3.0 
Bevel CF SISNIFICANCE = 0.05 


i ca EP meee ce ces ces we ee eee ee ee se Se eee See ee eee eee ee ea se 


N Magee! rere 2. ere 5 Pees | Shas 2 CON REG 
4 C.067 0. 050 0. 883 0. €42 C.580 Gg. 320 
€ Ce. O90 OQ. C93 0.817 0-689 0.719 0.945 
8 0.262 0.085 O.595 0.924 0.718 0.S56 

10 Ort 9 0.125 0.696 0.355 C.704 Cecat 

l2 Ga2o5 0.106 C.635S OP a C.&98 Conca 

14 Co 374 0.111 Gar Oag2c O./21 0.540 


> ce me ani i ce SE ec oe ee ee ee ee ee ee ee ee ee ee ee es es es ee ee ee eee 
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Bn Oo ae ee 


AXLTIT.VARTIA TION OF TYPES OF INTERVALS WITH WN 
Geek. © GEE . = O5 00 


SIGMA = 3.0 
maveleCr SIGNIFICANCE = 0.05 


= ee ce ee cee ee cee ee ee ae ee ee ee cee ee i i i i Mo i i td 


N moze! Tete? *. PyYeEes> “Roe sl (RSG 2 CONSE 
Cas O. C68 C. 824 0.713 0.750 Ooo 
0.226 0.105 Ol Aeie) 0.863 0.781 0.946 
0.376 O.118 0.506 Ones Geaecd Coco 
10 C2445 Geely 0.438 Oorc ao C2629 Gs 43 
Ne on iol Q- 110 Ceez g Oe clans 0.855 0.949 
14 0.57% OM Shee Oie2 25 Jae wee 0.857 0.548 
More VARIATION OF TYPES CF INTERVALS WITH % 
Gann. COEFF, = 0. /70C 
SIGMA = 3.0 
MEL SE SIGNIFICANCE = 0.05 
N Tyee: eferre 2° TYPE 3° PRGGe? Pees 2 CCN REC 
4 0.216 0.077 Ono 7 Os02G C.870 Cs oo2 
6 C354 0.122 0. 524 Jecos Caer, Ca oo 
8 Ce Dom Q.118 Caos 0.926 0.864% 0.943 
10 0.736 0.088 0.176 0.938 0.898 0.945 
iz 0-841 ORO GZ 0.097 0.941 dace 3s o $44 
14 0.3897 0. C24 C. 06S 0. 9590 Ce oiZ Cai SoZ 


<==. <a> <a ee cee eee ee es ee ee es ce es ee ee ee ee ee ee ee ee = oe ees ee es ee es ee es es ee ee ee ee ee ee es ee ee ee ee ee ee ee ee ee 
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XXVe VARIATION OF TYPES OF INTERVALS WITH AN 


COIR. CGEFF . = 9-300 
STGM4 = 5.0 
mevee GF SIGNIFICANCE = 0.05 


== ee ce we we se eee ee see > cmp eee CoP ee ce ae ee ee ee ee ee ee ee ee ee ee ee ee es ec ee ee ee ee ee ce 


N igewl. Uke 2. -fYPE 3S, Page) “eae o2 CON REG 
4 C.068 0.051 0.881 Genes 2 Gre tes 0.554 

6 Ce. 088 0. C94 C. 818 0.682 0.713 0.545 

8 0.138 0. 993 0.769 0.884 Q./10 ie Soh 
10 0.182 Operas | 0.691 0.363 Carat Cac3s 
2 Geo oS O. 107 0.628 Cao > C. 701 Cae 

14 C520 O. 110 C.540 0.920 Oeus 0.940 


= <oe See -e ee e e e eee ee  ee e eee e e  e e e  e e 


vice ven lati ON OF TYPES CF INTERVALS WITH N 


Caer COEFF . 
SIGMA = 39.0 
Weveu CE SIGNIFICANCE = 0.05 


\ 
© 
e 
mn 
© 
© 


= a cep oe ee eee ee eee es ee ee ee ee ee eee ee ee ce ee ee se ee ee ee ee 


N roe “TYPE 2 1YPE 93> eNes 1 ees 2 CON REG 
4 OSOYy O. C69 0.824 0.701 Om3 9 0.550 
OD 02109 0.676 0.360 GC. 7/61 ° S44 
C.363 0.125 0.511 Oe S2¢ C. 810 0. S5C 
10 Co 444 OP ae eS) C.441 Je FUL 0.817 OR ers) 
Vee 0.560 O.111 Os329 J 2953 0.856 0.949 
14 0.564 C2097 Oie259 0.950 0.845 Coe c 


<p es es es ee ee es es es se oe es es ee eee ee eee ee ee a > ca ee cee ee ee ee ee 
° 
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a Gs 2 #3’ ~~ - Oo a ~~. ee. hh) CLO 


XXVIIT.VARIATIOCN CF TYPES OF INTERVALS WITH N 


Coins c CEFF. = 0.900 
SIGMA = 5.0 
meee OF SIGNIFICANCE = 0205 


</> cee Oe ce ee ees ee eee ee es ee ce 2 <p oe oe ce ce es ee ee ee ee ee ee ee ee ee ee es ee ee re ee ee 


N mye 1” TYPE 2 TYPE 3 "PROB ') “FROB 2° CON REG 


> a cp ce ee eee ee we we a ee ae eee ee a ee eb ee EP ee eee oe eo 


10 0.901 0. 052 0.947 0.948 0.942 8510) 
eZ C.962 0.015 OP SON 2 0.947 C.938 C. S4E 
14 Cet oS 0.007 0.014 0.552 1.000 Speers. 


a cc A ec ee ee ee ee eee ee > oe ee ee eee ee Oe ee ee 


XXVITI.VARIA TION OF INTERVALS WITH SIGMA 


Beare CE SIZE 14 
eek .-COEFF. = Q. 7/00 
LEVEL CF SIGNIFICANCE = 0.05 


Ew a ee aE ce — 2 ee ee oe —_ ee Se a ee ee eee eee ee ee ee 


SIGMA Myeeel TYPE 2 Wee S PRB 1 PROS 2 CON SEG 


> a eee ee em eS ee ce eee ee eee ee Se 


= > a ee ee eee ee <a ce oe ee es ee eee eo a ee ee ee ee es es es ee es ee ee eee 
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XXIXe VARIATION GO INTERVALS WITH SIGMA 


SpaMePLE SIZE = 18 


CORR.CCEFF. = 0.500 
LEVEL JF SIGNIFICANCE = 0.05 


o_o ee ee ee ee ec eee ee ee ee ee ee ee ee ee ee es es eee ee es ee ees ee es es es es es es ee es ee es es ees ee ee ee ee ee eee ce 


SIGMA fictive coemen iy ress.) PROG Lee FeoB 2 CCN REG 
iD C2540 0.083 0.377 Oa Si2 C.819 €, coc 
2 Coeiers O- 116 0. 488 Oe 715 Ce Slo 0. S47 
3 C.376 O. 118 CE BNOls 0.923 0.322 O50 
4 0.372 0.118 0.510 0.922 0-822 0.$50 
> 0-363 0.126 0.511 0.928 C.810 CaS c 


a a es ce se se eee es se ee ee eee em eee a ce ee 


emcee boa do SE INTE®VatsS WITH SIGMA 


Serer Le BT Ze =22 
oa 2G JEFF « = 0.500 
Meet YF SIGNIETLCANCE = 0.05 


ee eee ee ee ee ee ee ee ee ee ee ee eee eee ee eee ee ee ee ee eee ee ee eee eee ee ee ee eee ee eee ee 


SIGMA Wace es 2 7 Y Meo PRIB Lt PRB 2 CON REG 
ih C.483 OQ. 110 0. 407 0. S05 Ceold C. S36 
2 Co 446 O. 115 C2439 0. 901 0.826 Gags 6 
3 0.445 O-11/ 0.438 0.899 0.828 GS 
+ 00446 0.115 0) ASR eee, Oa og? Ceesice Os 55 
5 0-444 0.115 0.441 0. 9U1 rerore ? OP ie 2. 


> cee ee ee ee es ee ee ee es es eee ee ee ees ees ee eee ee ee em ees es es es es eee ee es eee eee ee eee eee ee 


ay 





-_> 


PROB ORS iYVPESsOr SINT . 


LL, TYPES OF ENTERVALS VS p 





0 

6g 

8 

7 

6 

5 

4 

3 3, 

2 

: 2. ° 9 

My 5 3°45 6.8.9 1. 
CORR. COEFF. p > 

SAMPLE SIZE = 30 
SIGMA Sh eieag 
SIGNIFICANCE LEVEL = 0.05 


1, 2 and 3 on the curves indicate the types of Intervals. 
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> 


PROB. OF TYPES OF INT. 





2. TYPESRORe GER AeSe Vise 


iL 
0 
ju 
6 8 LO Wy igs 
N > 
CORR. SCORE? = 0.5 
SIGMA os MoO 
SIGNIFICANCE LEVEL = 0.05 


5) 









al 


PROB. OF TYPES OF INT. 


ce 


TYPES OF INTERVALS VS N 


6 8 LO gy 
N > 

GORE. CORFE. = 0,/7 

SIGMA =~ 1250 


SIGNIFICANCE LEVEL O05 


60 


14 





-> 


PROB. OF TYPES OF INT. 


4. TYPES OF INTERVALS VS o 





6 8 LO eZ 14 


CORR COPE er 0.) 


SAMPLE SIZE = 8 
SIGMA = 50 
SIGNIFICANCE LEVEL = O05 
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> 


PROB. OF TYPE 1 INT. 


oF 





PROD. On siya Eee Vo 
(VARIATION WITH a) 


= 0.02 
= Uns 
= 0, Le 
6 8 10 lls 14 
N + 
CORRE Onis = 0.9 
SIGMA = IAC 
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> 


Oa Oo 


LO 


30 


40 
50 


60 


10 


wine Vor CALCULATE PREBABILITY CF TYPE ly ™YPE 2sTY°E 
eet ceo ERVALS FOR STMULATED SAMPLES. 
PROGRAMMER T S WURTHY SEP 57 9.6 


ote a she sie oe oe oe ke ge oe ates te og teak ok ok oko ke kk KKK RAK AKKAE A AAR HSH 


ERMENS EGON Z2(55),X(55 ),V(100555),X1(50)+Y1(50),5(055), 


LIC(S)_IR( 5),19C(5),STAT(10, 7) sVS(0 109102 7) ,TVS(L0) 


CALL OVFLOW 

INC EX =1 

SIGMA=1.0 

ReEAD(S,2) Kot 

Weeire (6:2) K,T 
memeaAd( 1X,i2s 2X F 5-3) 
Wet kK seo. 0 } GO TO 460 
ROW=0.10 

Ba 309 13=1,5 

99 10 J=H1l95 
B=(1l-ROWs *2 ) ¥#.5 
SaeeeaArreN OF SAMPLES 
ISEED=12345 

IPC (J) =0 

CJ 250 [4=1,19 

03 5C M=1,100 

CALL SNORM(ISEED,Z;4«) 
X(L)=SIGMA*Z(1) 

C} 30 J=27K 

X(J )=20W= X( J~ 1)+3 =S IGMA*Z(J) 
O03 40 L=lyK 
VOM gL) =X01L) 

CIN TINUE 

0} 60 I1=1,5 

IR(I1) =0 

C2 z00 1=1,10C 

DQ 70 J=l,K 

S( J) =V¢(1,J) 

K-K=5 

KK=K/2 
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80 


$0 


100 


110 


J a1) 


DJ €C L=lyKK 

YI Cl) =8 (2 ¥*L) 

N=KK—-1 

DJ $0 LL=1,N 

Meeeev—t Ss (2b 1) +s (25b +1) )/2.0 
XSUN=0.0 

YSUM=0.0 

SXX=CG.0 

SXY=0 20 

SYY=0.0 

D9 100 KL=1,5 

IC (KL) =0 

99 110 M=1;N 
YSUM=YSUM+YTI(™) 
XSUM=XSUM 4X T() 

SONTINUE 

XB = XSUM/N 

YB8=YSUM/N 

D9 120 M=1,N 

SX X=S XX+( XI (M)—-X6) ** 

See Spier lM) =¥ 6) 2 

Sais oA tM )=X8) FOYT M)-Y¥8) 
CONTINUE 

VRES=(SYY -C(SXY¥*2 J/SXX)I/ (N-2) 
EH= SXY/SXX 

AH=Y8-BH* XB 

Soe 2 J*VRES 

Pee 2 I= SS/SXX 

O=S (27kK}—-YS 

== croh= P} 

C= (O%*2 )-(SS¥*(Nt1) 7M 

b= (ee teax 2 )*SXX/SS 

S50 = O—N-(P==2)/(SS*(Nt1)) 
Peete 5S 6) GM 79 150 
[er steel eO) GI) TC 500 
Daneeee XB) —S(K—-1) + 2. *BHEP/A 
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LS) 


ire © 
Lé0 
170 


200 


220 
250 


275 


2183) 
290 


E=( 2. /A)* (( ( SS*( PHX2) /SXX)+ (N41) ¥AXSS/N) *¥65 ) 
PIL=D-E 

PIR=D+E 

PVAL=S( K+ 1) 

Per ect. 1.0) 50 to 130 

IC (2 d= 1C(2)41 

IPOPVAL.LE-PIR .OR. PVAL.GE.PIL) 10(5)=1C(5)+1 
Gd TO 160 

IC (1 = IC(1) +1 

[EeeetERe VAAN «PIR. GEEp YA 16 G4) =1e(4)+) 
eeen 160 

fee n=10(3) +1 

Dy 1a) dese 

IR(J)=ER( J) 41 C(J) 

K=K+5 

CONTINUE 

me 220 Meee 

[PC ( J) SIPC J) #12 (J) 

CONTINUE 

ORINT STATISTICS 

STATCI9,1)=20h 

STAT(IB+2)=1PC{1) /1C00.6 

STAT (i8;3)=1PC(21/1000.0 

ST AT (18) 4 = 1PC(3 )/1000.0 
MEGoteCre; 2) .£0.0.C) Ga Ter27s 

STAT( IB 95)=1PC(4) /{ STAT(1B;2) #1000.0) 

GO TC 280 

ST AT (18,5 )=1°C(4) 

Ten r( 18,3) 260 ose) cGolta 2685 

STAT(IB 26) =1PC(5) /( STAT (18,3) *1000.0) 

Gemnc 290 

ST AT(1B6 J=19C(5 } 

SIG=STAT( 1B +2)#STAT( 18, 5)+STAT(I B73) *STATC IB; E) 


1 +STATCIB 94) 


STATCIB,7 J=S IC 
RVW=RIW+0 22 
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300 Ca pels 
DQ 305 19=1,5 
C9 305 J9=1,7 
VSCINDEX, 190 sJOI=STAT(IQ,JQ) 
20> CINTINUE 
IVS( INDEX )=N 
INDE X=INDEX+ 1 
IS Z=K-5 
WRITE(E,310) 182, SIGMA,N 
We EL O35: 325.) 
WRI TE( 69350) (CSTATOCK,L) plLaly7) »K=1_5) 
310 PORMWAICIXs! SAMPLE SIZE = ',15,' STIGMA = ',FS.0; ‘ 
IN = ', 15,7) 
eo FIRWAT( IX, * , OW TYPE L PY Gee 2 oiigs 
LFROCE.1L FRCB.2 CON. EG mes ee Oo te) 
350 FIRMAT( 7 F8 239 2X )y/ ) 
NR ITE( 6,485) 
a) RBS ll 
460 COR RK=0.1 
INCEX= INDEX—-1 
CJ 470 J=135 
WRITE (6,472) CORR, SIGMA 
WRITE(6 9475) 
DC 471 T=1,INDEX 
WRI TEC 6 e480) CIVS(L) pOVSCT 9d 9K) pK H257)) 
oat CIN @BNUE 
CIRF=ECIRR +4 22 
470 CIN TINUE 
Siz DORM ec OrescW@err, = ',F5.3,' SIGMA = ',F5.0;//) 
415 RPOeRment’ SAMPLE SIZE “YPEL Weyecce, Us ame, PQQ 
tee 1 ee BZ Cee aoe 4/7, 10l*=') } 
480 FIRMAT(I5,5X,€(F 8.392%) + /) 
500 STOF 
ENC 
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APL PROGRAM LINTREND 


V LTRTREMI 3X 
Sa REL Ch+O, 5X CU Ne RA) ) 
>?) MF EOS Chl 1 J 
C3] 9 Xe (e-2) 4X 
C4) KECFRYXOLS 
CS “alervd 
C6] Me] 
G7 J DPLLX UK CMeEM+] 
C8] MLE CXL ((2xM)—-3) +X (0 (62xXM)-1)7)xX0.5 
C9] MDL EMT gh} 
Oa, 7 
C119 rex C2xK KI 
Cea weKkK—1 
£13] Tlettari 


C143) [tHe VECTOR «I IS "3x2 
C15 Pte VECTOR Trl IS "3 rt 
C16] ' THE VALUE OF XF IS '3xF 


C173 Cierpy 

C18] C2e\H 

Ci9] C3ext 

C20] UMTE (Jet) PCL yC2,CZ 

C21] OMeRUuMT 

[227 SeUumMT +, x Om 

besd° SteR> 

G4 BHe (SI+,xUMT)+,xrt 

C25) SHS¢e(+/ (0 (C7TLT)«2))-CCORTI+, XOM) +, x BH) 
C26) AOeSTILCL Lathe xCC2xST CLs 27) + (h+1) «SIC 2527) 
bay deer O+1 
PesieeereCexetl 1s SI +2xStT Ley gaxCitt]) 
Pevyeeeoece tl Ss 33 

C30] CFe4xSHS+(h-3) 

eS eee Orta 5 x2) -CARKCE ) 

ener (ox ele Ss | XW CONE OX K I~ BWC 1 a - SRE 2a xK) )-ALXxcr 
C33] Ce (WxD)-AOXCF 

C34] 2IMx) (EF x2) L4xAxc 

C35] Fe ( (Bx2)-4xAxC)n0.5 

C36] &ek+A 

7 eee eer Ek C2xK p-j]j-eSA 


C38) ' THE LEFT LIMIT IS ';50'-E 
eS 4 ' THE RIGHT LIMIT IS ' 5046 
C40] 70 
C41] IM3' HMO INTERVAL EXISTS ! 

V 


oy) 





fd 
ie) 
Loe 
C4] 
ew 

Bo 
C7 I 

C8) 

Cae 

C altoua 
Cab aa 
a 
Coa 
C14] 
Cove 
BS] 
te? 
Css: 
eee 
e201 
C2008 
22 
peo 
C24] 
feo. 
E26] 
ie 7 J 
LAagele 
a? 
ESO J 
Cea ae 
L Bye 
Cae 
C34] 
sa J 
ESS 
7 2 
ESos 
ES 2 
C40] 
C41) 
C42) 
C43] 
C44] 
C45] 
C46] 


V 


Vv 


APL PROGRAM PREDICT 


Pee h oye to gtiels gi 9 (> ott ol gs SSR eS S36 s Stites OKT GS ay gp oSIe 5 Vs 
A oo IS ASSUME OO TO GE A VECTOR OF EVER HO OF ERMTE LES 
KECEHIXO.S 

mLe 1 O 

Me 1 

POKYK CMEME] 

per (2X Mie Soi OX Me) XO. 

MIEKLT y X41 

93 

TIEXLOx KY 

HeK—j 

TrertparyTr 

THe (+/ TI) =i 

MEE C/I) Hh 

SA Tée+E/ OC CKIAMKE) x CrI-rke) ) 

SMMeE+/ (CK IAKE) &D) 

STTE+/((TI-TE) 2&2) 

Khe Soy + SNK 

AH¢e TER~EH XE 

VRES¢e¢ (STT4+ O71 XO CSN 1D) HGH) ) PMH D 

Ses ZOGMaVOOS -SGeloe "24 //6 24357. 21447 ne. oOe 21300 
G¢eGsy, 2eeol LPH. Si 2,16 2145 ee eae ary ial 2 LOL 2.09 
GeG, 2,08 2.074 2.069 2.9064 2.906 2.056 2.052 2.048 2 
CeGy 2.041 24904 2+94 2+94 OO o8 RE LOSS IO COS 

SeCye Oe Ieee Old 2.02 2402, 2202 92,02 2,010 2:01 2401 
Te GLht-2] 

Aé (BHD) —( C=58%8%) xS¢( (TKD) XVRES ) ) 

Be 1X (DXBHXF Ee CHD OxXKJ-TE) ) 

CEC PRD) -SX( (MEL) EN) 

Fe (CBHRD) X SNS 

' THE F STATISTIC IS ee 

SSQt1—(Ma (M41) ) XCF a2) +5 

Jat 8 alee | 

SNR) OO ali 

3~T3X\F ¢SSQ 

SHeiterval IS OF Tike 2 

+) 
T14271-0 

De (OXKER HR (DKK) -]L J+ 2x BHXrF+A 

eer xe x Cok ( CCR ee SS ee Cie) it) tA) Ono 

' THE LEFT LIMIT IS 'S UE 

' THE FIGHT LIMIT IS 'SO+E 

' THE LENGTH OF INT IS calc xs 5 

+0 
Te toe fity 6 6IS 6 6OF TYPE 44 

+0 
Vet ee LT £S OF TYrFE 3 * 
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